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ABSTRACT 

The master equations describing processes of biological evolution in the framework of 
the random neighbor Bak-Sneppen model are studied. For the eqosystem of iV species 
they are solved exactly and asymptotical behavior of this solution for large N is analyzed. 
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1 Introduction 



The model of biological evolution proposed by Bak and Sneppen [1,2] describes 
mutation and natural selection of interacting species. It is the dynamical system 
that is defined as follows. The state of the ecosystem of iV species is characterized 
by a set {x\, xn} of N number, > Xi > 1. In so doing, X{ represents the barrier 
toward father evolution of the species. Initially, each Xi is set to a randomly chosen 
value. At each time step the barrier xi with minimal value and K — 1 other barriers 
are replaced by K new random numbers. In the random neighbor model (RNM), 
which will be considered in this paper the K — 1 replaced non-minimal barriers are 
chosen at random. 

The RNM is the simplest model describing the avalanche- like processes, which are 
supposed by a conception of " punctuated equilibrium" in biological evolution. These 
processes are the most characterizing feature for self-organized criticality recently 
intensively investigated both numerically and analytically [3-6] The RNM is more 
convenient for analytical studies. The master equations obtained in [3] for RNM are 
very useful for this aim. In [7] the explicit solution of master equations was found 
for infinite ecosystem. In this paper we solve the master equations for finite number 
N of species in ecosystem. We restrict ourselves with the simplest case K = 2. 

2 Master equations for RNM 

The master equations for the RNM are obtained in [3] . They are of the form: 

P n (t+1) (t) + B n+1 P n+l (t) + C n ^P n ^(t) + 

D n+2 P n+2 (t) + (£i<$ n , + A 1 S nA + Ci<$ ni2 )P (*) (1) 

Here, P n (t) is the probability that nat the time t is the number of barriers having 
values less than a fixed value A at the time t; < n < N, < A < N, , < t; P n (0) 
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are proposed to be given. For < n < N 

71 — 1 

A n = 2A(l-A) + ^- T A(3A-2), 

71 — 1 

5n=(l-A) 2 + ¥ - T (l-A)(3A-l), 

a = A 2 -^A 2 , D n ={l-Xf^- V (2) 

and A n = B n = C n = D n = for n = 0, n > N . 
In virtue of the definition of P n (t), 

Pn(t) > 0, (3) 

N 

T,Pn(t) = l. (4) 

n=0 

Making summation in (1) over n and taking into account (2) it is easy to establish 



that 



E^»(*+l) = E^n(*)- (5) 

n=0 n=0 

Therefore, if P n (0) are chosen in such a way that (4) is fulfilled for t — 0, then 
in virtue of (5) it is the case for the solution of (1) for t > too. For analysis of (1) 
it is convenient to introduce the generating function q(z,u): 

oo N 

q(z,u) = J2J: P n(t> n ^- (6) 

t=0 n=0 

In virtue of (3), (4) q(z,t) is polynomial in z, analytical in u for \u\ < 1 and 

= t^— (7) 
1 — u 

The master equations (1) can be rewritten for the generating function q(z,u) as 
follows: 

>•*> - = < J - * + A -~» 2 (K 1 - H (; " I)) x 

x(?(z,u)-g(0,«))+?(0,«)J. (8) 
The function 0) = J2n=o Pn(0)z n in (8) is assumed to be given. 



3 Asymptotic expansion of q(z, u) for big N 

If the function q(z, 0) has an asymptotic expansion in the region of big N of the 
form: 

then the equation (8) enables one to obtain the similar asymptotic expansion for 
q(z,u): 

q{z,u) = ^ 



to(N-l) k ' 

The main approximation of qo(z,u), the function q(z,u) can be found from the 
equation 

(z-u(l-X + Xz) 2 )q (z, u) = zq (z, 0) + u(l - X + Xzf{z - l)g (0, u) (9) 
following from (8). Since q (z,u) is analytical for \z\ < 1, \u\ < 1, 

= ag (Q;, 0) + u(l - A + Aa) 2 (a - l)g (0, u) (10) 

where 



Zi /A Hi 



is the solution of equation a — u(l — A + Aa) 2 = 0. Obviously, for sufficient little 
\u\, \a\ < 1. Thus, from equation (10) the function q o (0,u) can be found: 

/ n \ <?o(a,0) 

5o(0,u) = - • 12) 

1 — a 

Substituting (12) in the right hand side of (9), one can find the solution in the 
following form: 

„ ( v ..'i _ ^o( 2 ' + 1M1 ~ A + A^) 2 g (q, 0) 

go(2,wj — t j- : . , 2 , ( 7 , [i-o 

(z — u(l — A + Xz) z )(l — a) 

where a(u) is defined by (10). 
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For k > the functions qu{z, u) are defined by recurrent relations 

_ u(l - A + Az) 2 Qz - l)q k (g, 0) + (f - a)zq k (z, 0) 

(z-u{l-\ + \zf)(l-a) + 
t u(l — X + Az) 2 (r fe _i(z, u) + r fc _i (g,M)) 



z-w(l- A + Az) 2 
Here, 

9fi(z,u)-ft(0,)i) 

r k {z,u) = z- . (15) 

az z 

In virtue of (13), (14) the first correction to the lowest approximation (12) of q(z, u) 
has the form 

, _ u(l - A + \z) 2 {z - l)gi(a, 0) + (1 - a)zgi(z, 0) 
gi(Z ' Mj " (z- M (l-A + Az) 2 )(l-a) + 



■u(l — A + \z) 2 (r (z, u) + r (a, u)) 
+ A + Az) 2 ' 

where 



(16) 



3 (l-a)g (z,w) + (w(l-A + Az) 2 -l)g (a,0) 
ro( "' M) "^ (z- W (l-A + Az) 2 )(l-a) • (17) 



4 Exact form of q(z, u) 

Let us introduce the quantity 

Q(z,u) = . (18) 

z 

It follows from (8) that this function fulfills the relation of the form: 

(z — u(l — A + Xzf - -± ^— ^ " ^)<9(*, M ) = 

= g(z, 0) + (u(l - A + Az) 2 - l)g(0, u). (19) 

This inhomogeneous differential equation for Q(z,u) has a special solution 

i 

Q(z, u) = {N- l) e R(z ^ J e- R(x ' u) g(x, u)dx = Q sp (z, u). (20) 



(23) 



Here, 

R(z, u) = ^(ln(l - A + Xz) - (1 - u) ln(l - z) + ^-LzA—^ ), (21) 

_ g(x,0) + Kl-A + A^) 2 -l)g(0^) 

u(l-A + Az) 2 (l-z) ^ 
and the derivative of u) with respect to z has the form 

dR{z,u) _ (iV-l)(l-(l-A + A^) 2 ) 
dz ~ u(l- X + Xz) 2 (l- z) 

General solution of the corresponding to (19) homogeneous equation 

(, — u(l — A + Az) 2 - «d-A + A^(l- z) a )sfe u) = Q p4) 

is of the form 

Sfou) = F(u)e fl ^. (25) 

Here, F(u) is an arbitrary function of u. Hence, it follows from (19), (20), (25) that 
the function Q(z,t) can be represented as follows: 

Q(z, u) = F(u)e R ^ + Q sp (z, u). (26) 

In virtue of initial condition (7) for q(z, u) 

Q(l,«)=j4^-?(0,«) (27) 

For < u < l,z — > 1, S(z,u) diverges and S sp (z,u) has the finite limit: 

\imQ sp (z,u) = -i g(0,ii) (28) 

2— > 1 1 — M 

Hence, F(u) = in (25) and this representation for Q(z, u) can be rewritten in the 
form: 



Q{z,u) = {N-l)e R{z ' u) J e~ R(x ' u) g{x,u)dx + 

A 

A-l 
A 

+ (N - l)e R{z > u) J e- R ^g(x,u)dx (29) 



6 



It follows from (18) that 

n /A-l \(q(^,u)-q(0,u)) 
A ' ' A-l 

For the terms in the right hand side of (29) we have forO<u<l,0<A<l: 



lim e R{z ' u) = +00, (30) 

lim (N - l)e^»> [ e- R ^g(x, u)dx = X ^^[ u) ~ u) . (31) 
J A — 1 

Therefore (29) can represent the function Q(z,u) with necessary analytical proper- 
ties only if 

1 

J e -R(w) g ( XiU )dx = (32) 



A-l 
A 



This equation defines the function q(0, u): 



f p -R(x,u) q(x,0)dx 
A J X C (l-A+Ax) 2 (l-x) 

9(0, u) = (33) 

f r -R.(.r : n) (1-«(1-A+Aa) 2 )dx 
, J . (l-A+Ai;) 2 (l-x) 



A-l 
A 



Thus, we obtain from (29) the solution of equation (8) in the following form: 



A-l 
A 



A-l 
A 



where q(0,u) is defined by (33). 

5 Conclusion 



We constructed the solution of the master equation (8) for the finite number iV of 
species in the ecosystem. It can be proven that the main term (13) of its asymptotic 
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for large iV coincides with the one obtained in [7] . Using (34) one can obtain all the 
known analytical results for RNM. One can hope that it helps to understand better 
the most important properties of the self-organized criticality processes. 
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